Microbial signatures and enterotype clusters in fattening pigs: implications for nitrogen utilization efficiency.

As global demand for pork continues to rise, strategies to enhance nitrogen utilization efficiency (NUE) in pig farming have become vital for environmental sustainability. This study explored the relationship between the fecal microbiota, their metabolites, and NUE in crossbreed fattening pigs with a defined family structure. Pigs were kept under standardized conditions and fed in a two-phase feeding regime. In each phase, one fecal sample was collected from each pig. DNA was extracted from a total of 892 fecal samples and subjected to target amplicon sequencing. The results indicated an influence of sire, sampling period (SP), and sex on the fecal microbiota. Streptococcus emerged as a potential biomarker in comparing high and low NUE pigs in SP 1, suggesting a genetic predisposition to NUE regarding the fecal microbiota. All fecal samples were grouped into two enterotype-like clusters named cluster LACTO and cluster CSST. Pigs' affiliation with enterotype-like clusters altered over time and might be sex-dependent. The stable cluster CSST demonstrated the highest NUE despite containing pigs with lower performance characteristics such as average daily gain, dry matter intake, and daily nitrogen retention. This research contributes with valuable insights into the microbiome's role in NUE, paving the way for future strategies to enhance sustainable pig production.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits.

BACKGROUND: Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected. RESULTS: All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count. CONCLUSIONS: Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows.

BACKGROUND: Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS: In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION: Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Genomic analyses of nitrogen utilization efficiency, its indicator trait blood urea nitrogen and the relationship to classical growth performance and feed efficiency traits in a Landrace × Piétrain crossbred population.

Improving the nutrient efficiency in pork production is required to reduce the resource competition between human food and animal feed regarding diet components edible for humans and to minimize emissions relevant to climate or the environment. Thereby, protein utilization efficiency and its equivalent nitrogen utilization efficiency (NUE) play a major role. Breeding for more nitrogen (N) efficient pigs bears a promising strategy to improve such traits, however, directly phenotyping NUE based on N balance data is neither cost-efficient nor straightforward and not applicable for routine evaluations. Blood urea nitrogen (BUN) levels in the pig are suitable to predict the NUE and, therefore, might be an indicator trait for NUE because BUN is a relatively easy-to-measure trait. This study investigated the suitability of NUE as a selection trait in future breeding programs. The relationships to classical growth performance and feed efficiency traits were analysed as well as the relationship to BUN to infer the role of BUN as an indicator trait to improve NUE via breeding. The analyzes were based on a Landrace F1 cross population consisting of 502 individuals who descended from 20 Piétrain sires. All animals were genotyped for 48,525 SNPs. They were phenotyped in two different fattening phases, i.e., FP1 and FP2, during the experiment. Uni- and bivariate analyses were run to estimate variance components and to determine the genetic correlation between different traits or between the same trait measured at different time points. Moderate heritabilities were estimated for all traits, whereby the heritability for NUE was h2 = 0.293 in FP1 and h2 = 0.163 in FP2 and BUN had the by far highest heritability (h2 = 0.415 in FP1 and h2 = 0.460 in FP2). The significant genetic correlation between NUE and BUN showed the potential of BUN to be considered an indicator trait for NUE. This was particularly pronounced when NUE was measured in FP1 (genetic correlations r g = - 0.631 $$ {r}_g=-0.631 $$ and r g = - 0.688 $$ {r}_g=-0.688 $$ between NUE and BUN measured in FP1 and FP2, respectively). The genetic correlations of NUE and BUN with important production traits suggest selecting pigs with high growth rates and low BUN levels to breed more efficient pigs in future breeding programs.

Simulation of dual-purpose chicken breeding programs implementing gene editing.

BACKGROUND: In spite of being controversial and raising ethical concerns, the application of gene editing is more likely to be accepted when it contributes to improving animal welfare. One of the animal welfare and ethical issues in chicken breeding is chick culling, the killing of the male layer chicks after hatching due to the poor fattening performance. Although establishing dual-purpose chicken lines could solve this problem, unfavorable genetic correlations between egg and meat production traits hindered their competitiveness. Although it is also controversial in ethical terms, gene editing may accelerate genetic progress in dual-purpose chicken and alleviate the ethical concerns from chick culling. RESULTS: The simulation compared the utility improvement in dual-purpose use under two breeding schemes: one consisting in the improvement of the laying hens, and the second in the improvement of a synthetic line obtained from a layer broiler cross. In each breeding scheme, the breeding programs were simulated with and without gene editing. Polygenic breeding values and 500 simulated quantitative trait loci (QTL) with different levels of pleiotropy caused negative correlations between egg production, meat production, and overall health. The results of the simulation demonstrated that genetic gain could be accelerated by at most 81% for several generations if gene editing was used. The actual increase in genetic gain depended on the number of single nucleotide polymorphisms (SNPs) being edited per animal. The rate of genetic improvement became equal in scenarios with and without gene editing after 20 generations. This is because the remaining segregating QTL had small effects and their edition would have negative overall health effects from potential off-target edits. Although gene editing can improve genetic gain in quantitative traits, it can only be recommended as long as QTL with reasonable effect sizes are segregating and detectable. CONCLUSIONS: This simulation demonstrates the potential of gene editing to accelerate the simultaneous improvement of negatively correlated traits. When the risk of negative consequences from gene editing persists, the number of SNPs to be edited should be chosen carefully to obtain the optimal genetic gain.

Defining valid breeding goals for animal breeds.

BACKGROUND: The objective of any valid breeding program is to increase the suitability of a breed for its future purposes. The approach most often followed in animal breeding for optimizing breeding goals assumes that the sole desire of the owners is profit maximization. As this assumption is often violated, a generalized approach is needed that does not rely on this assumption. RESULTS: The generalized approach is based on the niche concept. The niche of a breed is a set of environments in which a small population of the breed would have a positive population growth rate. Its growth rate depends on demand from prospective consumers and supply from producers. The approach involves defining the niche that is envisaged for the breed and identifying the trait optima that maximize the breed's adaptation to its envisaged niche within the set of permissible breeding goals. The set of permissible breeding goals is the set of all potential breeding goals that are compatible with animal welfare and could be reached within the planning horizon of the breeding program. In general, the breed's adaptation depends on the satisfaction of the producers with the animals and on the satisfaction of the consumers with the products produced by the animals. When consumers buy live animals, then the breed needs to adapt to both the environments provided by the producers, and the environments provided by the consumers. The profit function is replaced by a more general adaptedness function that measures the breed's adaptation to its envisaged niche. CONCLUSIONS: The proposed approach coincides with the traditional approach if the producers have the sole desire to maximize their income, and if consumer preferences are well reflected by the product prices. If these assumptions are not met, then the traditional approach to breeding goal optimization is unlikely to result in a valid breeding goal. Using the example of companion breeds, this paper shows that the proposed approach has the potential to fill the gap.

The economic benefit of herd genotyping and using sexed semen for pure and beef-on-dairy breeding in dairy herds.

The cost benefits of herd genotyping and the benefits of using sexed semen have been affected by recent improvements in sexing technologies, the incorporation of direct health traits in the German total merit index for Holstein cattle, the deteriorating prices for purebred heifer calves and bull calves, and the introduction of herd genotyping programs. Inseminating the genetically superior dams with female-sexed Holstein semen increases the mean breeding value of the heifer calves and can produce more Holstein heifer calves than are needed for replacement. This provides an opportunity to increase the selection response in health and production traits at the farm level. A deterministic model is introduced that predicts the increase or decrease in net profit when a farmer takes part in a herd genotyping program and follows a certain insemination strategy. The types of semen that are allocated to cows and heifers are sexed semen or unsexed semen and Holstein semen or beef breed semen. The genetically superior heifers and cows are inseminated with female-sexed Holstein semen, intermediate dams with unsexed Holstein semen, and genetically inferior dams with unsexed or male-sexed beef breed semen. In general, participating in a herd genotyping program is beneficial for German Holstein breeders. The optimum proportions of cows and heifers that should be inseminated with a certain type of semen are sensitive to farm-specific peculiarities. A small price difference between crossbred bull calves and crossbred heifer calves often makes the use of male-sexed beef breed semen uneconomic. Under the conditions considered, it was found to be advantageous to inseminate approximately 50% of the heifers and 10% of the cows with the highest genetic merit with female-sexed Holstein semen. The optimum proportion of cows that should be inseminated with unsexed beef breed semen was found to be approximately 40%. In a herd with a low replacement rate, the selected heifers can exhibit their genetic superiority over a longer period of time, and a larger proportion of cows can be inseminated with beef breed semen. Participation in a herd genotyping program is, therefore, particularly beneficial for herds with low replacement rates.

Optimization Strategies to Adapt Sheep Breeding Programs to Pasture-Based Production Environments: A Simulation Study.

Strong differences between the selection (indoor fattening) and production environment (pasture fattening) are expected to reduce genetic gain due to possible genotype-by-environment interactions (G × E). To investigate how to adapt a sheep breeding program to a pasture-based production environment, different scenarios were simulated for the German Merino sheep population using the R package Modular Breeding Program Simulator (MoBPS). All relevant selection steps and a multivariate pedigree-based BLUP breeding value estimation were included. The reference scenario included progeny testing at stations to evaluate the fattening performance and carcass traits. It was compared to alternative scenarios varying in the progeny testing scheme for fattening traits (station and/or field). The total merit index (TMI) set pasture-based lamb fattening as a breeding goal, i.e., field fattening traits were weighted. Regarding the TMI, the scenario with progeny testing both in the field and on station led to a significant increase in genetic gain compared with the reference scenario. Regarding fattening traits, genetic gain was significantly increased in the alternative scenarios in which field progeny testing was performed. In the presence of G × E, the study showed that the selection environment should match the production environment (pasture) to avoid losses in genetic gain. As most breeding goals also contain traits not recordable in field testing, the combination of both field and station testing is required to maximize genetic gain.

Nitrogen and lysine utilization efficiencies, protein turnover, and blood urea concentrations in crossbred grower pigs at marginal dietary lysine concentration.

Nitrogen utilization efficiency (NUE) and lysine utilization efficiency (LUE) are key indicators of sustainable pork production and vary depending on nutritional and non-nutritional factors. The objective was to study NUE and LUE together with concentrations of blood urea nitrogen (BUN) and other metabolites in growing pigs fed diets with marginal Lys concentrations at 11-13 wk (40.5 kg mean BW) and 14 to 16 wk (60.2 kg mean BW). The cereal grain-soybean meal-based diets contained 10.6 and 7.9 g Lys/kg DM in periods 1 and 2, respectively. Feed intake and BW were measured for 508 individually penned pigs, and blood samples were collected 5 h after morning feeding at weeks 13 and 16. A subgroup of 48 barrows was used in a nitrogen (N) metabolism trial at weeks 13 and 16. In this subgroup, the mean N retention of pigs (27.3 g N/d) and mean LUE (70%) were not different between the periods, but NUE was higher in period 1 (47%) than in period 2 (43%) (P < 0.001). After administration of a single dose of 15N labeled glycine and measurement of 15N recovery in urine, the calculated whole-body protein turnover did not differ between the periods. The rate of protein synthesis was positively correlated with NUE (P < 0.001), but protein degradation was not. Excretion of urea-N in urine accounted for 80% of the total urinary N and was positively correlated with BUN. The N retention of all 508 pigs was estimated using an equation that was derived from the N metabolism data. N retention was on average 31.4 g/d, equal in both periods, and higher in barrows than in gilts in period 2, but not in period 1 (P = 0.003). The calculated NUE was, on average, 47% and was lower in barrows than in gilts (P < 0.001) and higher in period 1 than in period 2 (P < 0.001). The calculated LUE was, on average, 71%, and was lower in barrows than in gilts in period 2, but not in period 1 (P < 0.001). The BUN concentration was higher in barrows than in gilts (P < 0.001) and higher in period 1 than in period 2 (P < 0.001). BUN concentration was negatively correlated with NUE in Periods 1 (r = -0.50) and 2 (r = -0.15) (P < 0.05). We concluded that the maximum LUE was in the range of 70-72% under the conditions of this study, and only small differences between the periods and sexes existed. Protein synthesis, rather than degradation, appears to affect NUE. BUN concentration may be useful for estimating NUE in a large group of animals fed a diet with a marginal Lys concentration.

Increasing the utilization of protein and amino acids by pigs is an effective tool for improving the sustainability of pork production. This study showed that the utilization of lysine, the limiting amino acid in pigs, can reach 72% in growing pigs aged 11­16 wk when fed diets based on cereal grains and soybean meal. The average utilization of total nitrogen was 47%, with high variation among the pigs. Higher nitrogen utilization was associated with lower concentration of urea in the blood. This indicates that the blood urea concentration is an easy-to-determine proxy trait for nitrogen utilization in pigs when the concentration of lysine in the feed is low. Only small differences were observed between gilts and barrows during the investigation period. Nitrogen utilization increased when pigs had a higher rate of whole-body protein synthesis, which was associated with anabolic hormone concentrations in the blood.

Disturbed circadian rhythm of locomotor activity of pullets is related to feather pecking in laying hens.

Various aspects of activity, such as spontaneous activity, explorative activity, activity in open-field tests, and hyperactivity syndrome have been explored as causal factors of feather pecking in laying hens, with no clear results. In all previous studies, mean values of activity over different time intervals were used as criteria. Incidental observation of alternated oviposition time in lines selected for high (HFP) and low feather pecking (LFP), supported by a recent study which showed differentially expressed genes related to the circadian clock in the same lines, led to the hypothesis that feather pecking may be related to a disturbed diurnal activity rhythm. Hence activity recordings of a previous generation of these lines have been reanalyzed. Data sets of a total of 682 pullets of 3 subsequent hatches of HFP, LFP, and an unselected control line (CONTR) were used. Locomotor activity was recorded in pullets housed in groups of mixed lines in a deep litter pen on 7 consecutive 13-h light phases, using a radio-frequency identification antenna system. The number of approaches to the antenna system was recorded as a measure of locomotor activity and analyzed using a generalized linear mixed model including hatch, line, time of day and the interactions of hatch × time of day and line × time of day as fixed effects. Significant effects were found for time and the interaction line × time of day but not for line. All lines showed a bimodal pattern of diurnal activity. The peak activity of the HFP in the morning was lower than that of the LFP and CONTR. In the afternoon peak all lines differed with the highest mean in the LFP followed by CONTR and HFP. The present results provide support for the hypothesis that a disturbed circadian clock plays a role in the development of feather pecking.

Genomic dominance variance analysis of health and milk production traits in German Holstein cattle.

Genomic analyses commonly explore the additive genetic variance of traits. The non-additive variance, however, is usually small but often significant in dairy cattle. This study aimed at dissecting the genetic variance of eight health traits that recently entered the total merit index in Germany and the somatic cell score (SCS), as well as four milk production traits by analysing additive and dominance variance components. The heritabilities were low for all health traits (between 0.033 for mastitis and 0.099 for SCS), and moderate for the milk production traits (between 0.261 for milk energy yield and 0.351 for milk yield). For all traits, the contribution of dominance variance to the phenotypic variance was low, varying between 0.018 for ovarian cysts and 0.078 for milk yield. Inbreeding depression, inferred from the SNP-based observed homozygosity, was significant only for the milk production traits. The contribution of dominance variance to the genetic variance was larger for the health traits, ranging from 0.233 for ovarian cysts to 0.551 for mastitis, encouraging further studies that aim at discovering QTLs based on their additive and dominance effects.

Inferring causal structures of gut microbiota diversity and feed efficiency traits in poultry using Bayesian learning and genomic structural equation models.

Feed and phosphorus (P) efficiency are of increasing importance in poultry breeding. It has been shown recently that these efficiency traits are influenced by the gut microbiota composition of the birds. The efficiency traits and the gut microbiota composition are partly under control of the host genome. Thus, the gut microbiota composition can be seen as a mediator trait between the host genome and the efficiency traits. The present study used data from 749 individuals of a Japanese quail F2 cross. The birds were genotyped for 4k single-nucleotide polymorphism (SNP) and trait recorded for P utilization (PU) and P retention (PR), body weight gain (BWG), and feed per gain ratio (F:G). The gut microbiota composition was characterized by targeted amplicon sequencing. The alpha diversity was calculated as the Pielou's evenness index (J'). A stable Bayesian network was established using a Hill-Climbing learning algorithm. Pielou's evenness index was placed as the most upstream trait and BWG as the most downstream trait, with direct and indirect links via PR, PU, and F:G. The direct and indirect effects between J', PU, and PR were quantified with structural equation models (SEM), which revealed a causal link from J' to PU and from PU to PR. Quantitative trait loci (QTL) linkage mapping revealed three genome-wide significant QTL regions for these traits with in total 49 trait-associated SNP within the QTL regions. SEM association mapping separated the total SNP effect for a trait into a direct effect and indirect effects mediated by upstream traits. Although the indirect effects were in general small, they contributed to the total SNP effect in some cases. This enabled us to detect some shared genetic effects. The method applied allows for the detection of shared genetic architecture of quantitative traits and microbiota compositions.

Feed efficiency and phosphorus efficiency are of increasing importance in poultry breeding. It was frequently shown that next to the birds' genomes also the gut microbiota composition is important for these efficiency traits. The gut microbiota composition is a mediator between the genomes of the birds and their efficiency traits. In the present study, an approach was taken to consider the animal's gut microbiota diversity, efficiency traits, and the genomes of the animals together in a causal network to decipher the mediator role between the traits. Growing Japanese quail were used as model species. A stable network could be established that placed the diversity of the gut microbiota composition at the forefront, with direct and indirect links to other traits like phosphorus utilization and retention, feed per gain ratio, and growth. Together with genome scans, the results confirmed the mediator role of the gut microbiota composition because several traits associated variants affected the efficiency traits directly and indirectly via the gut microbiota composition.

Genomic regions underlying positive selection in local, Alpine cattle breeds.

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (FROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds.

Genetic and non-genetic factors influencing KLH binding natural antibodies and specific antibody response to Newcastle disease in Kenyan chicken populations.

This study aimed at investigating the influence of genetic and non-genetic factors on immune traits to inform on possibilities of genetic improvement of disease resistance traits in local chicken of Kenya. Immune traits such as natural and specific antibodies are considered suitable indicators of an individual's health status and consequently, used as indicator traits of disease resistance. In this study, natural antibodies binding to Keyhole Limpet Hemocyanin (KLH-NAbs) was used to measure general disease resistance. Specific antibodies binding to Newcastle disease virus (NDV-IgG) post vaccination was used to measure specific disease resistance. Titers of KLH-NAbs isotypes (KLH-IgM, KLH-IgG and KLH-IgA) and NDV-IgG were measured in 1,540 chickens of different ages ranging from 12 to 56 weeks. A general linear model was fitted to determine the effect of sex, generation, population type, phylogenetic cluster, line, genotype and age on the antibody traits. A multivariate animal mixed model was fitted to estimate heritability and genetic correlations among the antibody traits. The model constituted of non-genetic factors found to have a significant influence on the antibody traits as fixed effects, and animal and residual effects as random variables. Overall mean (±SE) concentration levels for KLH-IgM, KLH-IgG, KLH-IgA and NDV-IgG were 10.33 ± 0.04, 9.08 ± 0.02, 6.00 ± 0.02 and 10.12 ± 0.03, respectively. Sex, generation and age (linear covariate) significantly (p < 0.05) influenced variation across all the antibody traits. Genotype effects (p < 0.05) were present in all antibody traits, apart from KLH-IgA. Interaction between generation and line was significant (p < 0.05) in KLH-IgM and NDV-IgG while nesting phylogenetic cluster within population significantly (p < 0.05) influenced all antibody traits, apart from KLH-IgA. Heritability estimates for KLH-IgM, KLH-IgG, KLH-IgA and NDV-IgG were 0.28 ± 0.08, 0.14 ± 0.06, 0.07 ± 0.04 and 0.31 ± 0.06, respectively. There were positive genetic correlations (0.40-0.61) among the KLH-NAbs while negative genetic correlations (-0.26 to -0.98) were observed between the KLH-NAbs and NDV-IgG. Results from this study indicate that non-genetic effects due to biological and environmental factors influence natural and specific antibodies and should be accounted for to reduce bias and improve accuracy when evaluating the traits. Subsequently, the moderate heritability estimates in KLH-IgM and NDV-IgG suggest selection possibilities for genetic improvement of general and specific immunity, respectively, and consequently disease resistance. However, the negative correlations between KLH-NAbs and NDV-IgG indicate the need to consider a suitable approach that can optimally combine both traits in a multiple trait selection strategies.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle.

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Structural variation and eQTL analysis in two experimental populations of chickens divergently selected for feather-pecking behavior.

Feather pecking (FP) is a damaging nonaggressive behavior in laying hens with a heritable component. Its occurrence has been linked to the immune system, the circadian clock, and foraging behavior. Furthermore, dysregulation of miRNA biogenesis, disturbance of the gamma-aminobutyric acid (GABAergic) system, as well as neurodevelopmental deficiencies are currently under debate as factors influencing the propensity for FP behavior. Past studies, which focused on the dissection of the genetic factors involved in FP, relied on single nucleotide polymorphisms (SNPs) and short insertions and deletions < 50 bp (InDels). These variant classes only represent a certain fraction of the genetic variation of an organism. Hence, we reanalyzed whole-genome sequencing data from two experimental populations, which have been divergently selected for FP behavior for over more than 15 generations, performed variant calling for structural variants (SVs) as well as tandem repeats (TRs), and jointly analyzed the data with SNPs and InDels. Genotype imputation and subsequent genome-wide association studies, in combination with expression quantitative trait loci analysis, led to the discovery of multiple variants influencing the GABAergic system. These include a significantly associated TR downstream of the GABA receptor subunit beta-3 (GABRB3) gene, two microRNAs targeting several GABA receptor genes, and dystrophin (DMD), a direct regulator of GABA receptor clustering. Furthermore, we found the transcription factor ETV1 to be associated with the differential expression of 23 genes, which points toward a role of ETV1, together with SMAD4 and KLF14, in the disturbed neurodevelopment of high-feather pecking chickens.

Genomic rotational crossbreeding with advanced optimum contribution selection methods applied to simulated German Angler and German Holstein dairy cattle populations.

Many local dairy cattle breeds are facing genetic extinction due to a large proportion of foreign genes, which have been introgressed in the past. In addition, the performance gap to popular high-yielding breeds is increasing, resulting in a risk of numeric extinction. In the present simulation study, a genomic rotational crossbreeding scheme with the high-yielding German Holstein breed and the numerically small German Angler breed was analysed with the aim to utilize heterosis effects in the crossbred animals. Simultaneously inbreeding was controlled, and the amount of Holstein introgression observed in the Angler breed was reduced. Different scenarios of implementing OCS methods for Angler individuals were evaluated, which differed in their restrictions regarding kinship, native kinship, as well as the amount of genetic contributions from German Holstein. The results showed that rotational crossbreeding can result in superior crossbred offspring compared to the purebred parental lines, whereby OCS methods can simultaneously restrict the increase in inbreeding and keep the Holstein contributions at their current level. However, reducing the amount of migrant contributions while restricting the increase in the native kinship in Angler turned out to be a costly restriction. The reason was that Angler with low genetic contributions from Holsteins tended to have similar Angler ancestors. Consequently, reducing Holstein contributions would considerably increase the native kinship in Angler if it were not constrained. The constraint on the native kinship made a constraint on the conventional kinship superfluous and caused it to increase at a much lower rate than envisaged. This led to both, a high genetic diversity and a low genetic gain. The high genetic diversity in Angler also resulted in lower and oscillating heterosis effects in the crossbred animals. Thus, the reduction of migrant contribution did not increase heterosis effects in the crossbred offspring, and did not result in superior crossbred offspring in general.

Genetic analysis of production traits and body size measurements and their relationships with metabolic diseases in German Holstein cattle.

This study sheds light on the genetic complexity and interplay of production, body size, and metabolic health in dairy cattle. Phenotypes for body size-related traits from conformation classification (130,166 animals) and production (101,562 animals) of primiparous German Holstein cows were available. Additionally, 21,992, 16,641, and 7,096 animals were from herds with recordings of the metabolic diseases ketosis, displaced abomasum, and milk fever in first, second, and third lactation. Moreover, all animals were genotyped. Heritabilities of traits and genetic correlations between all traits were estimated and GWAS were performed. Heritability was between 0.240 and 0.333 for production and between 0.149 and 0.368 for body size traits. Metabolic diseases were lowly heritable, with estimates ranging from 0.011 to 0.029 in primiparous cows, from 0.008 to 0.031 in second lactation, and from 0.037 to 0.052 in third lactation. Production was found to have negative genetic correlations with body condition score (BCS; -0.279 to -0.343) and udder depth (-0.348 to -0.419). Positive correlations were observed for production and body depth (0.138-0.228), dairy character (DCH) (0.334-0.422), and stature (STAT) (0.084-0.158). In first parity cows, metabolic disease traits were unfavorably correlated with production, with genetic correlations varying from 0.111 to 0.224, implying that higher yielding cows have more metabolic problems. Genetic correlations of disease traits in second and third lactation with production in primiparous cows were low to moderate and in most cases unfavorable. While BCS was negatively correlated with metabolic diseases (-0.255 to -0.470), positive correlations were found between disease traits and DCH (0.269-0.469) as well as STAT (0.172-0.242). Thus, the results indicate that larger and sharper animals with low BCS are more susceptible to metabolic disorders. Genome-wide association studies revealed several significantly associated SNPs for production and conformation traits, confirming previous findings from literature. Moreover, for production and conformation traits, shared significant signals on Bos taurus autosome (BTA) 5 (88.36 Mb) and BTA 6 (86.40 to 87.27 Mb) were found, implying pleiotropy. Additionally, significant SNPs were observed for metabolic diseases on BTA 3, 10, 14, 17, and 26 in first lactation and on BTA 2, 6, 8, 17, and 23 in third lactation. Overall, this study provides important insights into the genetic basis and interrelations of relevant traits in today's Holstein cattle breeding programs, and findings may help to improve selection decisions.

Potential for quantifying general environmental resilience of dairy cattle in sub-Saharan Africa using deviations in milk yield.

Introduction: Genetic improvement of general resilience of dairy cattle is deemed as a part of the solution to low dairy productivity and poor cattle adaptability in sub-Saharan Africa (SSA). While indicators of general resilience have been proposed and evaluated in other regions, their applicability in SSA remains unexplored. This study sought to test the viability of utilizing log-transformed variance (LnVar), autocorrelation (rauto), and skewness (Skew) of deviations in milk yield as indicators of general resilience of dairy cows performing in the tropical environment of Kenya. Methods: Test-day milk yield records of 2,670 first-parity cows performing in three distinct agroecological zones of Kenya were used. To predict expected milk yield, quantile regression was used to model lactation curve for each cow. Subsequently, resilience indicators were defined based on actual and standardized deviations of observed milk yield from the expected milk yield. The genetic parameters of these indicators were estimated, and their associations with longevity and average test-day milk yield were examined. Results: All indicators were heritable except skewness of actual and standardized deviation. The log-transformed variance of actual (LnVar1) and standardized (LnVar2) deviations had the highest heritabilities of 0.19 ± 0.04 and 0.17 ± 0.04, respectively. Auto-correlation of actual (rauto1) and standardized (rauto2) deviations had heritabilities of 0.05 ± 0.03 and 0.07 ± 0.03, respectively. Weak to moderate genetic correlations were observed among resilience indicators. Both rauto and Skew indicators had negligible genetic correlations with both longevity and average test-day milk yield. LnVar1 and LnVar2 were genetically associated with better longevity (rg = -0.47 ± 0.26 and -0.49 ± 0.26, respectively). Whereas LnVar1 suggested that resilient animals produce lower average test-day milk yield, LnVar2 revealed a genetic association between resilience and higher average test-day milk yield. Discussion: Log transformed variance of deviations in milk yield holds a significant potential as a robust resilience indicator for dairy animals performing in SSA. Moreover, standardized as opposed to actual deviations should be employed in defining resilience indicators because the resultant indicator does not inaccurately infer that low-producing animals are inherently resilient. This study offers an opportunity for enhancing the productivity of dairy cattle performing in SSA through selective breeding for resilience to environmental stressors.

eQTL analysis of laying hens divergently selected for feather pecking identifies KLF14 as a potential key regulator for this behavioral disorder.

Feather pecking in chickens is a damaging behavior, seriously impacting animal welfare and leading to economic losses. Feather pecking is a complex trait, which is partly under genetic control. Different hypotheses have been proposed to explain the etiology of feather pecking and notably, several studies have identified similarities between feather pecking and human mental disorders such as obsessive-compulsive disorder and schizophrenia. This study uses transcriptomic and phenotypic data from 167 chickens to map expression quantitative trait loci and to identify regulatory genes with a significant effect on this behavioral disorder using an association weight matrix approach. From 70 of the analyzed differentially expressed genes, 11,790 genome wide significantly associated variants were detected, of which 23 showed multiple associations (≥15). These were located in proximity to a number of genes, which are transcription regulators involved in chromatin binding, nucleic acid metabolism, protein translation and putative regulatory RNAs. The association weight matrix identified 36 genes and the two transcription factors: SP6 (synonym: KLF14) and ENSGALG00000042129 (synonym: CHTOP) as the most significant, with an enrichment of KLF14 binding sites being detectable in 40 differentially expressed genes. This indicates that differential expression between animals showing high and low levels of feather pecking was significantly associated with a genetic variant in proximity to KLF14. This multiallelic variant was located 652 bp downstream of KLF14 and is a deletion of 1-3 bp. We propose that a deletion downstream of the transcription factor KLF14 has a negative impact on the level of T cells in the developing brain of high feather pecking chickens, which leads to developmental and behavioral abnormalities. The lack of CD4 T cells and gamma-Aminobutyric acid (GABA) receptors are important factors for the increased propensity of laying hens to perform feather pecking. As such, KLF14 is a clear candidate regulator for the expression of genes involved in the pathogenic development. By further elucidating the regulatory pathways involved in feather pecking we hope to take significant steps forward in explaining and understanding other mental disorders, not just in chickens.